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OPTIMIZATION  OF  ARMAMENT  DESIGN 


by 


H.  0.  Hartley,  R.  R.  Hocking, 
L.  R.  LaMotte  and  H.  H.  Oxspring 


1.  Introduction 

In  this  report,  we  describe  a  mathematical  programming  solution  to  the  problem 
of  optimizing  the  design  of  the  warhead  of  a  missile  for  the  purpose  of  destroying 
enemy  aircraft.  The  warhead  consists  of  fragments  and  explosive  charges.  Upon 
detonation  of  the  explosive,  the  fragments  are  projected  outward.  By  suitably 
packing  the  fragments  and  the  explosive  it  is  possible  to  control  the  fragment 
pattern,  thus  the  problem  is  to  determine  thst  fragment  pattern  which  will  be 
most  effective.  The  velocity  and  flight  path  of  the  target  relative  to  the 
intercept  missile  will  vary,  but  it  is  assumed  that  they  vary  according  to 
known  probability  laws.  In  view  of  this,  the  problem  to  be  solved  is  that  of 
determining  the  fragment  pattern  which  will  maximize  the  probability  of  destroying 
the  target. 

In  Section  2  we  develop  the  basic  notation  and  assumptions  for  the  special 
cose  in  which  the  flight  path  of  the  target  and  the  intercept  missile  are  in  the 
same  plane.  In  Section  3  the  optimization  problem  is  described  and  in  Section 
b  the  mathematical  programming  problem  is  formulated  for  this  coplanar  cose. 

The  formulas  necessary  for  implementing  the  solution  are  developed  in  Section  5 
and  an  example  is  given  in  Section  6.  The  generalization  to  the  ncn-coplonar 
ase  is  developed  in  Section  7,  and  other  extensions  developed  in  Section  8. 

2.  Notation 

The  notation  for  the  problem  will  be  developed  in  terms  of  Figure  1  which 
illustrates  t  typical  intercept  situation.  It  is  assumed  tha^  ot  time  t  =  0, 
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the  fissile  is  at  point  A  on  the  indicated  flight  path  traveling  with  velocity 

V,,  and  that  the  target  is  at  point  F  on  the  indicated  flight  path  and  has 

velocity  V^.  At  time  t  =  t  ,  the  missile  is  at  point  B  and  the  target  is  detected 

at  point  E.  The  angle  of  detection  0  is  known.  At  time  t  =  t^,  the  missile  is 

at  point  C  and  the  explosives  are  detonated  yielding  fragments  whose  velocity  is 

V  .  Finally,  a  fragment,  projected  at  an  angle  b,  is  assumed  to  intercept 
r 

the  target  at  point  D  with  t  =  t^.  The  angle  y  denotes  the  angle  at  which  the 

fragments  strike  the  missile  at  a  velocity  V  .  The  vulnerability  of  the  target 

s 

will  be  denoted  by  A(y,  V  )  which  is  assumed  to  be  a  known  function  of  y  and  V,. 

s  s 

The  following  quantities  are  assumed  known: 

(l)  VM,  the  missile  velocity 

(?)  t^  -  tp,  the  delay  time  after  detection 

(3)  Vf,  the  fragment  velocity 

(M  the  detection  angle 

The  velocity  of  the  target  V^,  the  angle  T)  between  the  two  flight  paths 

and  the  coordinates  of  the  point  E  are  not  assumed  known.  However,  it  is  assumed 

that  they  follow  a  known  probability  law.  This  probability  law  is  described 

in  terms  of  V_,  71  and  R  where  R  denotes  the  minimum  distance  between  the  mic,;i'’e 

r  mm 

and  the  target  if  they  follow  the  present  flight  path.  We  shall  denote  the  jci 
d  L"tr Ibution  of  VT,  T|,  Rm  by  f(VT>  71,  Rm) . 

?.  Foinulation  of  the  Optimization  Problem  in  the  Coplanar  Case 

The  probability  of  destroying  the  target  may  be  expressed  in  term-  of  R, 

the  distance  the  fragment  travel*;,  <?,  the  angle  at  which  the  fragment  is  pro,5  :t-d, 

the  angle  at  which  the  fragment  strikes  the  target  and  V  the  velocity  at 

s 

interception.  For  given  values  of  R,  P,  y  and  V  ,  this  probability  is  assumed 

s 


to  be  given  by 
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P(R,  9,  Y.  V8)  -  1  -  exp{-  A(y ,  Vs)  p(6)/  2}  .  (l) 

R 

The  function  p(6)  represents  the  fragment  density  function  in  terms  of  the  angle 

6.  (It  is  assumed  that  the  fragment  pattern  is  symmetric  about  the  axis  of  the 

missile.)  The  'intercept  geometry',  (R,  6,  y,  v  )  depends  on  the  missile  constants 

s 

V„,  t  -  t  ,  V  and  0  and  on  the  target  variables  V_,  Th  and  R  which  are  dis- 
M  a  L)  r  T  m 

tributed  according  to  f(V„,  71,  R  ).  In  Section  5,  the  formulas  for  determining 

T  m 

the  intercept  geometry  (R,  6,  y,  V  )  in  terms  of  the  missile  constants  and  the 
target  variables  are  developed.  Thus  a  probability  distribution  is  induced  on 
the  intercept  geometry.  Denote  the  joint  distribution  of  the  intercept  geometry 
variables  by  g(R,  6,  y,  V  ).  The  average  probability  of  destroying  the  target 
is  then  given  by 

Jg(R,  9,  V,  V6)  {l  -  exp(  -  A(y,  Vs)  p(6)/R2)}  (2) 

where  the  integration  is  over  the  range  on  the  variables  (R,  6,  y,  V  ). 

s 

The  objective  is  to  determine  the  density  function  p(0)  so  as  to  maximize 
the  'kill'  probability  given  by  (2).  Since 

Jg(R,  Y.  Vg)  -  1  (3) 

it  is  sufficient  to  consider  the  minimization  of 

Jg(R,  6,  Y »  Vs)  exp  (  -  A(y»  Vg)  p(6)/R?)  .  (M 

The  function  p(®)  is  constrained  to  satisfy  the  requirements 

Jp(  Q)d9  =  M  (5) 

and 

(O 


p(6)  >  0 


for  all  6. 
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The  constraint  (5)  is  the  requirement  that  a  given  number  of  particles  will  be 
packed  into  the  warhead  and  (6)  is  the  obvious  requirement  that  the  density 
function  be  non-negativ^. 

The  minimization  of  (4)  with  respect  to  p(6)  where  p(6)  must  satisfy 
(5)  and  '6)  is  a  problem  in  control  theory,  that  is,  the  calculus  of  variations 
further  complicated  by  the  inequality  constraint  (6).  Only  special  cases  of 
this  general  problem  have  been  solved  and  the  complexity  of  the  integrand  in  (4) 
forces  us  to  use  an  approximate  method  of  solution. 

4.  The  Mathematical  Programming  Solution 

In  order  to  solve  the  problem  formulated  in  Section  3  it  is  necessary  to 
assume  that 

(i)  the  distribution  f(VT>  T),  R^)  is  discrete  and, 

(ii)  it  is  sufficient  to  determine  p(6)  for  only  a  finite  set  of  values  ar.l 
then  approximate  p{0)  by  a  LaGrange  interpolating  polynomial. 

With  these  assumptions  the  problem  described  by  (4)  (5)  and  (6)  is  now  written 

as 

minimize  E  gt  exp{  -  A(yt,  V^)  P^}  (7) 

subject  to  E  a  p  =  M  (8) 

i  1 

Pt  >  0  (9) 

The  sum  in  (7)  and  (8)  is  over  all  points  (R,  0,  y,  V  )  which  hove  positive 

s 

probability  in  the  discrete  distribution  g(R,  9,  y,  V  ).  The  variables  p. 

S  1 

represent  tne  value  of  the  fragment  density  function  at  the  corresponding  values 
of  0,  that  is,  p^  =  p(0i).  The  constraint  (8)  is  just  a  suitable  eppreximation 
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to  the  constraint  (5)  for  the  relevant  values  of  6. 

Rather  than  determine  the  density  at  the  points  0^  as  determined  by  the 

distribution  g(R,  6,  y,  V  )  it  has  been  decided  to  specify  in  advance  a  set  of 

s 

variables  8*,  k  =  1,  r  at  which  p(0)  will  be  determined.  The  approximation 

to  p(0)  will  then  be  given  by  a  LaGronge  interpolation  over  these  points,  that  is 


where 


p(£)  =  I  LR(e)pk»  Pk  -  p(9£) 
k=l 


V6>  ■ 


n  (0  -  e*) 

q=l  q 

_ 

S  (6*  .  e*) 

q=i 

qA 


(10) 


(ID 


The  functional  (7)  is  then  written  as 


E  G 
1 


i 


(12) 


and  the  minimization  is  over  the  pv. 

A 

It  now  remains  to  determine  the  distribution  g(R,  0,  y,  Vg)  induced  by  the 

distribution  on  V^,  T),  R^.  The  appropriote  formulas  will  be  developed  in  Section 

b.  In  fact,  it  is  not  necessary  to  explicitly  develop  this  distribution  but 

rather,  the  functional  (12)  is  expressed  directly  in  terms  of  the  variables 

V_,  T),  R  os  follows: 
i  m 


N 

E  f 
J*1 


i 


exp| 


A(V  vs 


V  V"j}  • 


(13) 


In  (13),  it  is  assumed  that  the  distribution  f(VT,  T),  R^)  assigns  probability 


to  N  points,  (V  ,  71,  R  )  of  f  ,  j  =  1,  N.  Further,  R  ,  b  ,  y  ,  V 

1  m  J  J  J  J  5j 

denote  the  values  cf  these  variables  for  a  given  point  (V^,  T|,  ^m)j  as  computed 
from  the  formulas  in  the  next  section. 

The  function  (13)  is  seen  to  be  convex  in  the  variables  and  hence  its 
minimization  subject  to  the  constraints  (8)  and  (9)  is  just  a  problem  in  convex 
programming  and  can  be  solved  by  a  number  of  algorithms.  In  Section  6  we  describ 
an  example  using  the  Hartley-Hocking  (  _  .^)  ecr.v.  ::  programing  algorithm. 


Development  of  Formulas  for  the  Coplanar  Case 


The  solution  of  the  programming  problem  developed  in  Section  U  required  that 

\:e  develop  expressions  for  R,  0,  y  and  Vg  in  terms  of  V^,  7),  Rm  as  well  as  the 

constants  V  ,  t  -  t  ,  V_  and  ft. 

M  a  U  F 

For  convenience,  we  shall  assume  that  t^  ■  0,  hence  the  point  (Xq,  Yo)  re  ore 

cents  the  coordinates  at  which  detection  is  made.  In  general,  the  point  on  the 

target  which  we  wish  to  hit,  say  the  target  center,  will  not  be  the  same  as  the 

point  which  is  first  detected.  Let  (X*,  Y*)  be  the  coordinates  of  the  target 

'•enter  at  t  =  0  where  X*  «  X  +  c_ ,  Y*  =  Y  +  c.. 

o  o  1’  o  o  2 

Let  the  coordinates  of  the  target  center  at  time  t  be  denoted  by  (X^.Y^)  whr 


xt  =  x*  -  t  VT  cos  7) 
Yt  *  Y;  -  1  vt  sln  ”  • 


Similarly  the  coordinates  of  the  missile  (X^,  Y^)  at  any  time  t  are  given  by 

v  v 


<-tVM 


The  coordinates  of  a  fragment  projected  at  angle  6  for  any  time  t  >  t^  are  given 
by  (<•  <>  where 

<  *  ‘  v„  *  -  V  VF  9 

<-(t-  td)  VF  sin  e  (1G) 


0  <  e  <  2n  . 
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The  first  task  in  describing  the  intercept  geometry  is  to  determine  the  time 
of  intercept  t^.  This  is  achieved  by  solving  the  equations 


Using  (14)  and  (16)  the  equations  (17)  become 


X*  -  t  V.  cos  T)  -  t  VM  =  (t  -  t ,)  cos  0 

o  t  '  M  d  F 

(18) 

Y*  -  t  VT  sinT]  =  (t  -  td)  Vp  sin  6  . 

Squaring  both  sides,  adding  and  simplifying  yields  the  following  quadratic  in  t: 

t2(V?  +  4  *  2  VTVM  COS  11  - 
-2t(XJVT  cos  I)  *  x;vM  ♦  YJVt  sin  n  -  v£td) 

+  <X1  *  *  0  d9) 


For  certain  values  of  the  input  variables  and  missile  constants,  this  quad¬ 
ratic  may  have  either  zerc,one  or  two  reel  roots  which  exceed  t^.  The  usual  case 
will  be  that  only  one  root  exceeds  t^,  and  this  will  be  t^.  If  no  roots  exceed 
tH  for  the  case  that  is  (VT,  T],  R  )j  then  it  is  impossible  to  intercept 


.th 


with  this  choice  of  input  variables.  In  this  case,  the  j  term  is  deleted 

from  (13)  and  the  solution  obtained  os  usual.  However  the  kill  probability 

(2)  will  now  be  decreased  by  an  amount  f  .  The  case  of  two  roots  greater  than 

J 

t^  indicates  that  two  'hits'  are  possible.  At  the  present  time  we  only  consider 
one  and  hence  let  t^  be  the  smallest  root  greater  than  t^. 

The  angle  ft  yielding  the  interception  is  ntw  determined  from  (l8)  by  setting 
t  **  tj  and  solving  for  6.  The  distance  the  fragment  must  travel,  Rj,  is  then 


given  by 
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Er  ■  <xtx  -  v/ +  • 

Finally,  the  angle  at  which  the  fragment  strikes  the  target  is  given  by 


where 


t8n  y*  = 


V  =  V_  sin  9  +  V_  sin  T1 

sy  F  T  1 

V  =  V„  +  V_  cos  9  +  Vm  cos  71  . 

sx  M  F  T  1 


The  striking  velocity  is  given  by 


=  V2  +  V2 
x  sy 


Inspection  of  these  re  alts  shows  that  we  have  expressed  the  intercept 

geometry  (R,  8,  y>  V  )  in  terms  of  X*,  Y*  the  coordinates  of  the  target  center  at 

s  o  o 

t  =  0  rather  than  the  miss  distance  R  as  desired.  The  modification  of  the 

m 

formulas  is  easily  achieved.  For  a  specified  detection  angle  0,  the  detection 
point  satisfies  the  equation 

Yq  =  XQ  tan  0  .  (23) 

The  squared  distance  between  the  missile  and  the  detection  point  at  any  time  t  is 


given  by 


R2  =  (Xo  -  t(VT  cos  T)  +  VM))2  +  (Y0  -  t  VT  sin  7))2  .  (24) 


This  distance  is  seen  to  be  minimum  for  t  =  t  where 

m 


X  (V_  coo  T|  +  V.)  +  Y  V_  sin  T\ 
t  _  ov  T _ 1  M  o  T  1 

m '  i  *  it +  2  Vm  cos  ” 


(25) 
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In  view  of  (23)  we  have 


t 


m 


VT  cos  71  +  VM  +  VT  sin  7]  tan  0 

VT  *  VM  ♦  ?  VtVH  C0S  ’I 


=  Xo  H(Vt,  Vm,  71,  0)  . 


(26) 

(26) 


The  square  of  the  miss  distance  is  thus  given  by  substituting  (23)  and  (26) 
into  (2h).  Thus, 

Rm  =  Xo  I*-1  “  ^VT  C0S  ^  +  VM>H]2  +  ^tan  ^  "  VT  Sin  71  H-*2} 

*  G(Vt,  Vm,  7),  0)  (27) 


Thus,  for  a  specified  value  of  R  we  can  determine  X  and  hence  X*  and  thus 

m  o  o 

specify  the  intercept  geometry.  The  procedure  is  summarized  as  follows: 

(i)  From  (27),  determine 

X  =  R  <f*  and  X*  *  Xn  +  c.  (28) 

o  m  o  o  1 

(Note:  It  is  assumed  that  G  is  different  from  zero.  The  case  G  =  0 

VT  sin  7] 

arises  if  tan  0  =  ;; — ■  and  corresponds  to  zero  miss 

VT  cos  71  +  VM 

distance. ) 

(ii)  From  (23)  determine  Yq  and  Y*  =  Yq  +  c0  and  then  determine  the  time 
of  intercept  t^  from  (19). 

(iii)  The  intercept  geometry  (R,  0,  y,  V  )  is  then  determined  from  (20), 

s 

(18),  (21)  and  (22). 

Thus  each  term  in  the  functional  (13)  may  be  evaluated.  Using  the  some 
information  the  derivatives  of  (13)  with  respect  to  p  msy  be  evaluated  for  any 

lv 

particular  point  =  p*.  These  two  evaluations  are  all  thst  are  necessary  in 


11 


the  Hartley-Hocking  convex  programing  algorithm  and  hence  the  optimal  solution 
is  readily  obtained. 


6.  Illustration  of  Computations  in  Coplanar  Case 


The  optimization  problem  as  formulated  in  Section  I4  required  the  minimization 
with  respect  to  p^f  k  =  1,  r  of 


J/j  'xp{-A<v 


(29) 


subject  to  the  constraints 


E  atpk  =  M 
k=l  *  k 


(30) 


p  >  0  k  =  1,  r 


(31) 


The  p^,  k  =  1,  r  represent  the  value  of  p( 6)  for  specified  angles  9*, 

k  =  1,  r  and  the  approximation  to  p(9)  is  then  given  by 


P(°)  -  ^Ve)pk 


(32) 


where  L^(  6)  are  the  LeGrange  interpolation  coefficients  as  given  by  (ll).  It 
is  clear  that  the  constraint  (31)  is  not  sufficient  to  guarantee  that  the 
approximate  fragment  density  function  (32)  will  be  non-negative  for  all  8. 

To  partially  rectify  this  situation  we  adjoin  to  the  problem  the  constraints 


E  ^ Vpk  - 0  J  s  •••> N  • 


(33) 


Since  these  constraints  are  linear  in  the  variables 


pk 


the  problem  may  still  be 
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solved  using  the  convex  programming  algorithm.  In  addition  to  giving  added 
assurance  that  the  approximation  to  p(  6)  will  be  non-negative,  the  constraints 
(33)  insure  that  the  optimal  solution  of  (29)  is  between  zero  end  one  which  is 
essential. 

As  observed  earlier,  the  minimization  of  the  functional  (29)  subject  to 

the  constraints  (30) ,  (31)  and  (32)  by  the  Hartley-Hocking  algorithm  requires 

that  we  be  able  to  evaluate  (29)  and  its  derivatives  with  respect  to  p  for 

k 

specified  values  of  p  ,  k  =  1,  r.  In  this  section  we  shall  illustrate  the 

A 

procedure  used  to  prepare  the  problem  for  solution.  In  particular,  consider  the 
evaluation  of  one  term  of  the  functional  (29),  that  is 

fj  exp{  *  A(v  vSj>  V9.tVRj}  •  (3M 

Recall  that  the  missile  constants  Vu,  t ,  -  t_,  V_  and  0  are  fixed  but  that  N 

M  a  Dr 

triples  (VT,  T|,  Rm)jf  J  =1,  • . . ,  N  are  specified. 

The  example  considered  has, 

V„  =  5000  ft. /sec. 

Vp  =  11,000  ft. /sec. 

A(Y.,  V  )  «  1000 

3  j 

t .  -  t_  =  .002  sec. 
d  D 

0  ■  n/u 

the  density,  p(6),  will  be  determined  for  the  points 
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The  distributions  on  V^,  T), 


and  R  are  assumed  to  be 
m 


r;vT)  .  { 
f(1)  »  { 

f<V  ■  { 


i 

0 

1 

3 

0 

1 

3 

0 


VT  =  UOOO,  U250,  U500,  5000  ft. /sec. 

otherwise 

n  „  n 

^  =  °»  "  5 

otherwise 

Rq  =  50,  100,  150  ft. 
otherwise. 


The  numbers  used  for  this  test  problem  are  not  bssed  on  real  data  but  are  merely 
designed  to  provide  a  simple  example. 

For  simplicity,  it  is  assumed  that  the  variables  V^,  T)  and  R^  are  independent, 

thus 

f(VT,  1),  Ra)  3  f(VT)  f (T|)  f(Rn)  . 

The  first  step  in  the  evaluation  of  r2'^  is  to  determine  which  triple  (V  ,  T],  R  ) 

1  HI 

is  being  considered.  Let  V_  =  1*500,  T|  3  0,  R  =  100.  We  s^r  that  there  r.r- 

T  ra 

M  =  2r  triples  and  ir.  this  sirr.plo  example, 

f  -  ^  J  =  1,  ....  N  .  (35) 

We  then  determine  Xq  according  to  (28)  and  Yq  according  to  (23)  as 

X  =  100 

0 

Y  =  100  . 

0 

V’e  assume  X  =  X*  and  Y  =  Y*  for  this  example.  I4-  should  he  observed  that,  ir. 
o  o  c  o 

general,  detection  could  occur  for  either  f  =  n/,  <-r  "f  =  — n/,  .  w«  adopt  th« 

*4  4 

convention  that  tar.  /  >  0,  -,'r.et  is,  y  >  0. 


-  lU  - 


Substituting  into  the  quadratic  (19)  and  solving  yields 

tj  =  0.0111 

We  now  solve  (18)  for  0.  Let 


sin 


e  = 


Y 


o 


-  VT  tT  sin  7) 

Ih  -  VVF 


v 


cos 


(VT  COS  7)  +  VM)tr 

-  ‘A 


The  value  of  0  is  obtained  according  to  the  following  rules: 

(a)  If  arc  sin  (v)  <  0  then 

(1)  0  =  2tt  +  arc  sin  (v)  if  cos  9  >  0 

(2)  0  »  tt  -  arc  sin  (v)  if  cos  9  <  0 

(b)  If  arc  sin  (v)  >  0 

(1)  0  =  arc  sin  (v)  if  cos  0  >  0 

(2)  9  =  tt  -  arc  sin  (v)  if  cos  0  <  0. 

In  this  case  we  obtain  the  intercept  angle  corresponding  to  the  current  choice 
(VT,  71,  Rm)j  as 

0j  =  1.626  RAD  . 

From  (lU)  the  position  of  the  target  at  intercept  is  given  by 


Xtr  *  Xo  -  tIVT  COS  1 

=  100  -  (45CO) (0.0111)  =  50.5 

Y.  =  Y*  -  tV,  sin  7)  =  100 

t_  o  T  1 
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and  hence  from  (20)  we  obtain 


2  k 

R.  =  1.1602  x  10  . 

3 


Next,  we  must  evaluate  the  LaGrange  interpolation  coefficients  1^(6)  for  the 


current  value  of  8  =  1.6257.  For  example, 
J 


9 

n  (0.  -  e*) 

_  _p  j 

Ll(9j)  =  9 - -  °*oco'3 

V6!  -  V 


Repeating  this  for  k  =  2,  3,  . r  =  9  would  ./held  the  terms  1,^(6^), 
k  =  1,  ...,  9- 

Finally,  from  (21)  and  (22)  we  could  obtain  y.  and  V  and  evaluate 

0  j 

A(y,  V  ) .  In  the  current  exainple  the  function  A(y,  V  )  was  assumed  constant 
s  s 

co  this  is  not  necessary.  At  this  point  we  have  available  the  necessary 
information  to  evaluate  (3*0  for  any  choice  of  p.  and  also  its  derivative 
with  respect  to  any  p  .  The  procedure  described  above  is  then  repeated  for 

A 

all  possible  choices  (V^,  T),  to  generate  the  N  terms  in  the  sum  (29). 

It  should  be  emphasized  that  all  of  this  computation  is  done  prior  to  the 
application  of  the  convex  programming  algorithm.  This  algorithm  will  not  be 
described  at  this  time. 

The  approximation  (30)  to  the  integral  constraint  (5)  was  trapezoidal. 

Thus 


a.  =  1  k  =  2,  . . . ,  8  . 

k  ’ 


The  value  of  M  wes  assumed  to  be  2000  — . 
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For 
in  Table 


this  data  the  convex  programming  algorithm  gave  the  solution  shown 

1. 


TABLE  1.  Solution  to  Example  Problem 


k 

pk 

1 

0.250  TT 

0 

2 

0.3125  n 

0 

3 

0.3750  TT 

0 

4 

0.4375  n 

1359.2 

5 

0.5000  n 

433.9 

6 

0.5625  tt 

17.8 

7 

0.6250  n 

109.9 

8 

0.6875  tt 

79.2 

9 

0.7500  n 

0.0 

The  approximation  to  p(8)  is  then  given  by  (10)  and  (ll).  The  kill 
probability  (2)  is  given  by 

P  =  1  -  0.1535  =  0.8U65 

Execution  time  on  the  IBM  360-65  for  this  example  was  0.0^  minutes. 
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7.  Development  of  the  Three-Dimensional  Case 

In  this  section  we  consider  the  extension  of  the  problem  to  the  more  general, 
three-dimensional  case.  A  typical  intercept  situation  is  illustrated  in  Figure 
2.  The  X-oxis  is  again  the  flight  path  of  the  intercept  missile.  At  time 
t  =  0  the  target  is  detected  at  position  (XQ,  Yo>  Zq)  traveling  with  velocity 
VT  on  n  flight  path  described  by  the  angles  T)  and  to .  At  time  t  =  =  0  the 

target  is  detected  with  detection  angle  0.  At  time  t  =  t^  the  explosives  are 
detonated  and  finally  at  time  t  =  t  a  fragment  projected  in  a  direction  described 
by  the  angles  a  and  R  is  assumed  to  intercept  the  target.  The  striking  angle  y 
and  velocity  Vs  are  defined  as  in  the  planar  case. 

Following  the  development  in  Section  5  we  shall  now  develop  the  expressions 
for  the  variables  describing  the  intercept  geometry  in  terms  of  the  missile 
constants  and  the  target  variables.  Let  the  coordinates  of  the  target  at  time 
t  be 


where 


X,  =  X*  -  t  V_  cos  Ci  cos  T) 
t  o  T 

Y,  =  Y*  -  t  V_  cos  b  sin  T] 

O  T 

Z  =  Z*  -  t  V„,  sin  0- 
to  T 

X*  =  XQ  +  cx,  Y£  =  yQ  +  cp,  Z*  =  Z D  +  c^  . 


The  coordinates  of  the  missile  are 


‘  ‘  VM 
Vt  '  0 


(35) 


(36) 
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The  coordinates  of  a  typical  fragment  are  given  by  (t  >  t.), 


xf  *  t  V  ♦  (t  -  td)  Vpcos  B  cos  or 
vf=(t  -  t  )  V  cos  9  sin  o 

l  u  r 

z[  =  (t  -  td)  VF  sin  3 


0<o<2n  _^<0<| 


As  before,  the  time  of  intercept  tj.  is  determined  by  setting 
F 

Z  =  Z  ,  and  solving  for  t. 

L  X# 

Using  (35)  end  (37)  we  obtain 


_  V  - 

*t>  Xt  "  V 


XS  -  t  VT  cos  a;  cos  T1  =  t  VM  +  (t  -  td)  VF  cos  0  cos  a 

Y*  -  t  V_  cos  a  sin  7]  ■  (t  -  t  )  V_  cos  0  sin  a  (38) 

01  a  r 

Z*  -  t  VT  sin  0  =  (t  -  td)  Vp  sin  0 
Eliminating  o  and  0  yields  the  following  quadratic  in  t 

t?(VT  +  VM  +  2  VMVT  C0S  C0S  71  ’  VP 

-2t(X*VT  cos  0:  cos  7]  +  X*VM  +  YJVt  cob  U.’ Bin  71  ♦  Z*VT  sin  0  -  Vptd) 

2  2  2  ?2  .  . 

♦ xs  ♦ y;  ♦  K  -  vf  ■ 0  (39> 

The  last  real  root  of  (39)  which  exceeds  t  is  the  time  of  intercept  t^.  We  then 
obtain  Qj  and  0^  from  (38)  and  finally  is  obtained  from 

BI  *  \  -  v/  ♦  Tj  ♦  zlz  •  <'■-> 

Inspection  of  these  results  shows  that  we  have  described  the  intercept 
geometry  in  terms  of  X*,  Y*,  Z*  and  the  two  angles  o  and  0.  Assuming  that  the 


fragment  pattern  is  symmetric  with  respect  to  the  axis  of  the  missile,  it  is  not 
necessary  to  specify  O'  and  6  but  only  the  angle  that  the  fragment,  direction  make 
with  the  missile  axis.  Denoting  this  angle  by  6  as  in  the  planar  case  we  see 
that  0  is  given  by 


cos  6  =  cos  a  cos  3. 


(41) 


The  direction  of  the  fragment  and  its  velocity  relative  to  the  target  are 
determined  by  the  quantities 

V  -  V..  +  V„  cos  9  cos  O'  +  Vm  cos  <■’  cos  T) 
sx  M  F  T  1 


V  =  V„  cos  0  sin  a  +  v  cos  a'  sin  T) 
sy  F  T  1 


(4?) 


Vgz  =  Vp  sin  0  +  VT  sin  (•'  . 


The  angles  and  y?  as  indicated  in  Figure  2  are  given  by 

Yx  =  Vj  "  1)»  V2  *  Y|  - 


where 


V  V 

tan  y*  =  and  tan  y*  =  t— ^ 

Vsx  *  ” 


(v2  +v?)i 

V  sx  sy/ 


(**3) 


(44) 


The  squared  velocity  is  given  by 


■2  ”2  +  v2  +  V2 


V  =  V 


sx 


(>*5) 


We  now  consider  the  problem  of  expressing  the  intercept  geometry  in  terms  of 
the  minimum  miss  distance. 

For  a  specified  detection  angle  0  we  see  that 


P  ?  P  P 

Y  +  Z  =  tan  0  X 

00  o 


(*•6) 


where  (X^,  Y  ,  Z^)  is  the  detection  point.  The  squared  distance  between  the 


ion  point  an]  the  missile  rjt  any  time  t  is  given  by 
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R2  =  (XQ  -  t(VM  +  VT  cos  o-  cos  T|))2 

+  (Yq  -  V  t  cos  O'  sinT])2  (47) 

+  (ZQ  -  t  VT  sin  a;)2 


This  distance  is  seen  to  be  minimum  for  t  =  t  where 

m 


t  = 
m 


xo(vM  +  Vf  cos  o'  cos  'll)  +  Yo(Vt  cos  O'  sin  Tl)+ZQVT  sin  it' 

m 


VM  +  VT  +  2  VMVT  C0S  C°S  11 


The  square  of  the  minimum  miss  distance  is  thus  given  by  substituting  (48)  into 

(47).  To  use  the  formulas  developed  thus  far  we  now  attempt  to  determine 

X  ,  Y  ,  Z  in  terms  of  R  .  It  is  seen  that  sn  additional  angle  must  be  specified 
o’  o’  o  m 

namely  the  angle  which  describes  deviation  from  the  coplanar  case.  Accordingly, 
we  define  the  angle  ♦  by 

Z 

tan  =  y"  •  (^9) 

o 


Using  (49)  we  obtain 


Y  =  X  tan  0  cos  ♦ 
o  o 

(50) 

Z  =  X  tan  0  sin  ♦ 
o  o 

Substituting  (50)  into  (48)  yields 

»m  -  *0  h(vt.  vh,  -n.  U,  <f)  (51) 


and  hence  (47)  and  (51)  yields 


_  22  _ 

Rn  =  Xo  ‘  H<vm  +  vt  cos  U'  cos 

2 

+  (tan  f  cos  ^  -  VT  H  cos  0.' sin  7|)  (5?) 

+  (can  0  sin  ^  -  H  VT  sin  a.)^  . 

It  should  be  observed  that  there  is  no  loss  of  generality  in  assuming  that 

the  angle  ♦  is  zero  or  equivalently  that  Zq  is  zero.  This  is  equivalent  to 

assuming  that  the  X-Y  plane  in  Figure  2  is  determined  by  the  line  of  flight 

of  the  intercept  missile  and  the  point  at  which  the  target  is  detected.  Care 

must  be  taken  so  that  the  angles  7),  and  u  defining  the  target  flight  path  are 

described  relative  to  this  convention. 

In  view  of  the  above  comments,  we  see  that  either  X  ,  R  or  for  that  matter 

’  o  m 

the  variable  R  defined  as  the  distance  between  missile  and  target  at  the  time  o 
o 

detection  may  be  used  to  determine  the  intercept  geometry  keeping  in  mind  that 

all  of  these  variables  refer  to  the  point  of  detection  of  the  target.  The  choir 

of  which  of  these  variables  to  use  will  depend  on  the  data  available. 

From  this  point,  the  procedure  is  just  as  in  the  planar  case.  That,  is, 

specification  of  the  missile  constants  end  the  variables  VT,  7),  OC,  and  either 

R  ,  X  or  R  enables  us  to  determine  the  intercept  geometry  (R,  0,  y,  V  )  and 
moo  o 

hence  the  mathematical  programming  algorithm  may  be  used  to  determine  the 


optimum  fragment  density. 


6.  Extension  to  Multiple  Targets 

In  Sections  5  and  7,  the  formulas  for  computing  the  intercept  geometry 
for  the  special  cose  of  a  single  point  target  were  developed.  These  formulas 
were  then  used  to  develop  the  mathematical  programming  problem  for  the 
determination  of  the  optimum  fragment  pattern.  Frequently,  the  incoming  air¬ 
craft  will  have  more  than  one  vulnerable  point  or  target,  for  example  the 
engine  and  the  pilot  are  normally  thought  of  as  distinct  point  targets.  Usually, 
destruction  of  some  combination  of  these  point  targets  is  sufficient  to  destroy 
the  aircraft.  Such  combinations  as  (i)  two  out  of  three  (ii)  at  least  one  and 
(iii)  either  targets  1  and  2  or  targets  3  end  4  are  encountered.  The  development 
of  the  appropriate  mathematical  program  for  solving  this  extended  problem 
usually  causes  no  difficulty  but,  in  some  cases,  the  resulting  functional  may 
no  longer  possess  the  desired  convexity. 

To  indicate  the  procedure,  consider  the  special  case  of  an  aircraft  with 
two  point  targets  such  that  the  destruction  of  either  will  result  in  destruction 
of  the  aircraft.  Assuming  that  the  location  of  the  targets  (l  and  2)  relative 
to  the  point  of  detection  is  known,  then  using  the  formulas  of  Section  5  (or  7), 
two  intercept  geometries  (R,  9,  y,  V  ),  and  (R,  9,  y,  Vo)0  may  be  computed.  The 
associated  kill  probabilities,  say  and  are  given  by  (l).  The  probability 
of  destroying  at  least  one  of  the  targets  for  this  intercept  geometry  is  given  by 

1  -  (1  -  Px)(l  -  P2)  (53) 

Assuming  that  A,(y,  V^)  and  AD(y,  V  )  represent  the  vulnerability  functions  for 

X  3  d  S 


the  two  targets  then  the  function  to  be  minimized,  that  is,  one  minus  the 
average  hill  probability  (analogous  to  (13))  is 


-  2k  - 


N  r 

expl-^Yy,  Vg  )  ^(^j) 


4-1  J 


lj  k=l 


J* 

•  «P{-A3(Y2j>  %  > 


(5’0 


This  function  is  seen  to  be  convex  in  the  variables  p  and  hence  the  minimi- 

it 

zation  of  it  subject  to  the  constraints  (8)  and  (9)  as  well  as  (33)  is  a 
convex  programming  problem. 

More  general  requirements  on  the  combination  of  point  targets  which  must  be 
destroyed  to  insure  the  destruction  of  the  aircraft  cause  no  difficulty  in  formu¬ 
lation  but  in  general  lead  to  a  lack  of  uniqueness  and  optimaltiy  of  the  solution. 
For  example,  suppose  that  the  destruction  of  both  targets  is  required.  The 
probability  of  destroying  the  aircraft  for  a  given  intercept  geometry  is  given 
by  P,,  where  and  are  given  by  (l)  in  terms  of  (R,  0,  y,  Vp)^  and 
(R,  0,  Y,  V  )Q.  In  this  case,  after  reformulation  as  a  mathematical  program 
we  see  that  the  function  to  be  minimized  is  no  longer  convex.  In  this  cose, 
more  that  one  relative  minimum  may  exist.  The  algorithm  of  Fiacco  and 
McCormick  (19&M  may  be  used  to  determine  relative  minima  and  by  using  several 
initial  points  more  than  one  relative  minimi  m  ay  -V. •  y i i.  fho*  is, 
designer  may  have  several  fragment  patterns  suggested  each  of  which  are  locally 
optimal.  Assuming  that  they  ere  all  physically  and  economically  practical  he 
would  then  select  the  one  which  yields  the  highest  average  kill  prolability. 

There  is  no  guarantee,  however,  that  there  does  not  exist  a  superior  solution. 


References 


1.  Fiacco,  A.  V.  and  G.  P.  McCormick  (196*0  "Computational  Algorithms  for 

the  Sequential  Unconstrained  Minimization  Technique  for  Nonlinear 
Programming,"  Management  Science,  X,  No.  2,  601-617. 

2.  Hartley,  H.  0.  and  R.  R.  Hocking  (1963)  "Convex  Programming  by  Tangential 

Approximation,"  Management  Science,  IX,  No.  k,  600-612. 


Vctiritv  CUssifict'tvn 


DOCUMENT  CONTROL  DATA  •  R  &  0 

(Security  cts»  si  first  o(  MUt,  My  of  mkHrtef  snd  snnotMtlon  m  u»t  W  (WUrt^  whtn  th*  mrmrwll  ft 

V  OHiCciA  r.NO  activity  fC«CW^«  1#«  RIRORT  tk'CUHi  Ty 

Unclagolfled 

Tax#*  A&M  University  <».  aeove 

Unclassified 


MM  «tr*ral I  taper!  It  elattlllt.'l  | 

T  tSCCURlTV  CI.A5tl*l-:»TlOH  [ 

_ _ 1 C1.J  V 


»  rin.x 


OPTIMIZATION  OP  ABttAMENT  DESIGN 


*■  OCICMlMTi**  noth  a/  rtpart  mM  Inalvaiaa  MiImJ 

Final  Report 


»  WTnomiii  name,  mIMIa  Initial,  laai  turn*) 

H.  0.  Hartley,  R.  R.  Hocking, 

L.  R.  LaMotre  and  H.  H.  Oxaprlng 


July  1969 


M.  CON  TRACT  OR  GP 


yp.  total  no.  or  radii  |?t.  mo.  or  RIPS 


•  •a  ORIGINATOR**  AIPORT  NUM9LRIS) 


N00014-68-A-0140 

t.  MNOJtCT  NO. 

NR047-700 


Number  13 


tt.  OTHKR  PIC  PORT  MOHI  (Any  Qtht  t Udtnk9t§  After  mmy  •«  1/4*104 
Aft/#  rmpmri) 


*•.  MITRIIUTlOM 


This  docuarmt  has  bar.-  approved  for  public  release  and  sale; 
lta  distribution  La  unlimited. 


tfl.  IPDNSORINC  MI4.ll 


Office  of  Naval  Research 


In  this  report  a  method  la  developed  for  determining  the  optimal  warhead 
design  for  an  Intarcept  micelle.  Iha  problem  poecd  Is  that  of  determining  the 
fagreant  density  pattern  which  will  maximise  the  average*  probability  of 
destroying  an  incoming  aircraft  when  the  flight  parameters  of  the  tergne  are 
not  assvned  known  but  are  assumed  to  follow  a  known  probability  lnw.  The 
problem  is  formulated  aa  a  mathematical  pio.^rraxinp  problem  requiring  hi 
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